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Abstract. We study quantum transport properties of an open Heisenberg XXZ 
I ^~ spin 1/2 chain driven by a pair of Lindblad jump operators satisfying a global 

'microcanonical' constraint, i.e. conserving the total magnetization. We will show that 

this system has an additional discrete symmetry which is particular to the Liouvillean 

1 ^ ' description of the problem. Such symmetry reduces the dynamics even more than what 

would be expected in the standard Hilbert space formalism and establishes existence of 

' multiple steady states. Interestingly, numerical simulations of the XXZ model suggest 

£h that a pair of distinct non-equilibrium steady states becomes indistinguishable in the 

2 thermodynamic limit, and exhibit sub-diffusive spin transport in the easy- axis regime 

qH of anisotropy A > 1. 
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1. Introduction 

Simulating interacting quantum many-body systems in-, and in particular, out-of- 
equilibrium is at the forefront of the current experimental and theoretical research (see 
e.g. PQ, [2])- While it has been recognized that treating large isolated quantum-many 
body systems (with generic physical properties) poses a formidable task to any method 
of theoretical description or simulation, it has soon become apparent - first within the 
community of quantum optics [3] - that often a macroscopic number of observables can 
be treated as quantum noise, and that the behavior of a few essential physical observables 
can be extracted in terms of the formalism of open quantum systems [U [5]. From 
mathematical-physics point of view most notable is perhaps the discovery of Lindblad, 
Gorini, Kossakowski and Sudarshan [HI U\, that within the Markovian approximation 
any time-development of an open system's (reduced) density operator pit), acting on 
a N— dimensional Hilbert space % describing the (slow) degrees of freedom of interest, 
can be described by the quantum Liouville equation in the form 

N 2 -l 

p{t) = Cp(t) ee -i[H,p(t)] + £ (L m p{t)Ll - -{LlL m , p(t)}) (1) 

where [A, B] = AB - BA, {A, B} = AB + BA, H is the Hermitian Hamiltonian 
operator, and L m is a set of at most iV 2 — 1 so-called Lindblad operators, all over "H. 
We use units in which h = 1. C can be understood as a (super)operator acting on the 
space B(7i) of bounded operators over Hilbert space 7i. The space B(T-L) shall also be 
considered as a Hilbert space itself, equipped with the Hilbert-Schmidt inner product. 
The Liouvillian flow (fTl) is the most general one which is (i) local-in-time, (ii) respects 
the positivity pit) > 0, and (iii) conserves the trace tr p{t) = 1, and forms the so-called 
quantum dynamical semi-group. However, the approach of open-quantum-system has 
only quite recently been applied to the non-equilibrium problems of condensed matter 
physics [HI El EH [HI H21 EH CHI US] , such as the controversial quantum transport problem 
in one-dimension [TBI [T7J [18] . 

The central object of concern in physics applications is the fixed point of dynamical 
semigroup, p^ = lim^oo p(t), or the null- vector of Liouvillian 

tpoo = (2) 

and which shall be termed non- equilibrium steady state (NESS). This state is relevant 
in the asymptotic long time limit of the evolution of the system, and does not decay 
despite the dissipative nature of the Lindblad equation. The interesting case, when the 
fixed point is a pure state p^ = \ipoo)(4>oo\ is of particular appeal in quantum optics 
and quantum information theory [19], where such a state l^oo) is called a dark state. 
Sufficient condition for existence of a dark state is that |?/>oo) is an eigenstate of H and is 
annihilated by all jump operators, L^il)^) = 0. Dark states are also known as examples 
of decoherence-free states, and are of importance in quantum computing [20] . 

We note that the Lindblad equation (jlj can in general describe the system 
coupled to several thermal (or chemical, magnetic) baths at different values of the 
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thermodynamic potentials, hence the resulting asymptotic states p^ can be considered 
as intrinsically nonequilibrium. 

Of central importance is the understanding and control of uniqueness of NESS, 
for example for applications in quantum memories [21] which could be realized in 
terms of possibly degenerate (non-unique) NESSs. Another potentially very iteresting 
application of non-unique (multiple) steady states of Liouvillean quantum dynamics 
could be in detecting non-equilibrium quantum phase transitions, or regions of criticality 
|24j . Importantly, multiple fixed points of dynamical semigroups are also a crucial 
part of the concept of decoherence-free subspaces [20], one of the key ideas for fighting 
errors in quantum computation. Precise conditions for the uniqueness of NESS, i.e. 
independence of p^ of the initial condition p(0), has been established a while ago 
[22] [23] • Essentially, it follows from the Evans theorem [22J that NESS is unique if 
and only if the set of operators {H, L m , m = 1, 2, . . .} generates the entire multiplicative 
operator algebra B(7i). One the other hand, the situations where NESS is not unique 
have been systematically much less explored. For example, in the case of existence of a 
discrete symmetry, of either (i) Liouvillian £ as a whole, or (ii) of all members of the set 
{H, L m , m — 1, 2, . . .} one might expect non-unique NESSs classified by the eigenvalues 
of the symmetry operation. 

In this paper we will analyze how the existence of symmetries can allow us to 
treat the transport problem in a boundary driven open anisotropic Heisenberg XXZ 
spin 1/2 chain with a micro-canonical constraint, i.e. enforcing the exact conservation 
of total spin (magnetization). Quite interestingly, our numerical results suggest that 
the two NESSs corresponding to the zero magnetization sector characterized by the 
eigenvalue of the parity-like symmetry become indistinguishable in the thermodynamic 
limit, and that in the easy-axis regime (of anisotropy A > 1) the spin transport becomes 
sub-diffusive (insulating in the thermodynamic limit), quite different from the results 
|25| 126] for the un-constrained boundary-driven XXZ chain. This seems consistent with 
a recently claimed co-existence of diffusive and insulating transport in a gapped XXZ 
model [21] • I n the Appendix we provide a general mathematical prescription how in 
both cases (i,ii) of the previous paragraph the symmetry can be facilitated to reduce 
the problem and block-diagonalize the matrix of the Liouvillian superoperator, whereas 
only in the second case (ii) the non-uniqueness of NESS is guaranteed by the existence 
of a symmetry operator and distinct Liouvillian fixed points can be labelled by the 
eigenvalues of the symmetry operation. 

2. Symmetric boundary driven open XXZ spin chains 

In order to study a simple system which does not have a unique NESS and to provide 
some additional insight into the debate on quantum transport in one dimension [161 I2Z] , 
we shall discuss a finite open anisotropic Heisenberg XXZ spin 1/2 chain on n sites, 
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with the Hamiltonian 

n-\ 

H = J2 [°t°Zn + <°U + ^M + x] ■ (3) 

i=l 

o~* ,y,z designate a set of standard Pauli matrices corresponding to physical sites i = 
1,2, ... ,n and acting on a Hilbert space % = (C 2 )® 71 of dimension iV = 2™. We couple 
the chain to magnetic reservoirs only through the boundary spins i = 1 and i = n. 
Two cases can be considered: 

(i) Firstly, we take a set of four Lindblad operators 

lj = v/f(i - JT)*t, lj = Vra + ^Kt, (4) 

Ll = y/T(l + ix)<rl , LJ = y/T(l-fj)a- . 

where two real parameters T > and \x G [0, 1] designate, respectively, the strength 
of coupling to a pair of magnetic reservoirs and strength of (non-equilibrium) 
driving, af = (of ± iof )/2 are the spin-flip operators. 

(ii) Secondly, we consider a pair of Lindblad operators 

L\ = r(l - fi)afa~, LI = T(l + u)^a+, (5) 

again parametrized by two coupling parameters T, \x having the same interpretation 
as above. We note that the Liouvillian flow ^ equipped with ^ strictly conserves 
the total magnetization M = X/ILi a li a ^ t ne expense of non-local coupling between 
the first and the last site. In fact, we can imagine Lindblad operators here as 
incoherent quantum jumps which transfer spin excitations between i — 1 and 
i = n sites while conserving total number of spin-excitations (quasi-particles). 
Alternatively, such a model can be interpreted as a spin ring where hopping is 
fully coherent on bonds (1, 2), (2, 3) . . . , (n — 1, n) while it is fully dissipative (and 
asymmetric for \i ^ 0) on one bond (n, 1), and is a particular case of quantum 
exclusion process complementing the one presented in Ref. [28] (see also [29| 130]). 

We note that the quantum transport models described by the Hamiltonian pi) and 
Q, or (p)]), can be considered as out of equilibrium steady state models with minimal 
incoherent input, namely with incoherent processes taking place only on the boundary, 
or on the single bond. 

Let P be a permutation operator which exchanges the site i with n — i + 1 for all 
i. In fact, P can be written as a unitary operator over % with explicit representation 
in the computational basis (eigenbasis of a\, of |ai, 02 • • • a n ) = (1 — 2aj)|ai, 02 • • • a n ) for 
Oj G {0, 1}) as 

P= ^2 \ai,a 2 ,...,a n )(a n ,a n - 1 ,...,a 1 \ (6) 

(ai,a 2 ,...,a„)e{0,l}™ 

and can be interpreted as a parity operator. Combining it with spin flips on all sites we 
define another unitary (parity-like) operator 5" as 



)i. 



s = pY[°f- ( ? ) 



»=i 
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Clearly, S is a Z2 symmetry of the Hamiltonian, as [H, S] = and S 2 = 1. 

In case (i) with Q one notes that SL^tf = SL^ = LJ and SL% = LJ hence 
the Liouvillian flow (fl) exhibits a weaA; symmetry in the sense of Appendix, Eq. (11). 



Precisely the same case of symmetrically boundary driven XXZ chain (3|4) has been 



discussed in several recent papers [3TJ E2J ESI E31 EH ES] and has been shown to exhibit 
interesting non-equilibrium trasport properties. In this case, one can use Evans theorem 
[22J in order to prove that NESS is unique. The number of invariant subspaces is 
ns = 2, with symmetry S eigenvalues s\ = +1 and S2 = — 1. Let the complete set of 4 
symmetry adapted operator sub-spaces be now suggestively denoted as {s a , sp} = B a> p. 
As pointed out in the Appendix (also explaining the notation), the symmetry S can still 
be used in order to reduce the Liouvillian to two blocks E\ = {+1,+1}©{ — 1, — 1}, and 
£>2 = {+1,-1} © {— 1, +1}, which have roughly similar dimensions for large n. 

Let us now focus on case (ii). S now becomes a strong symmetry (Appendix, Eq. 



(12)) of the flow (1|5) as SL s m = L s m for m = 1,2. In this case, however, we have an 



additional strong (continuous, U(l)) symmetry, generated by magnetization operator 



Sm = e 1</,A/ , as also SmL^ = L s m and SmH = H. We shall say that the flow (l|5) obeys 



a micro- canonical constraint. Thus we can use the Theorem 4.1 of Appendix in order to 
prove the existence of a pair of distinct NESSs for each eigenvalue of total magnetization 
s z G {— n, —n + 2, . . . , n — 2, n}, i.e. we have a 2N + 1 dimensional convex set of fixed 
points of the Liouvillian flowj|] As from the quantum transport point of view the most 
interesting is the zero magnetization sector, we will in the following fix s z = and 
assume the size n to be even. There we have a line of degenerate NESS parametrized 
by a number from a unit interval 116 [0, 1] 

poo{u) = upl + (1 - «)/4, pi e {+1, +1}, pi e {-1,-1}. (8) 

2.1. Numerical results and transport properties 

Let us now turn to the physical properties of our microcanonically constrained open 
XXZ chain, in particular to the steady state spin current and magnetization profiles. 
The spin current operator corresponding to the bond (i,i + 1) is defined as, 

ji = afaj +1 - ajaf +1 = 2i(a;a+ +1 - a+a m ) (9) 



-i) 



and obeys the operator continuity equation of the form ^(<r|/2) = i[H, <r z /2] = ji —j 
for i = 2,3. . . ,n — 1. Steady state expectation values shall be denoted as J = (ji) 
(which does not depend on site index i due to continuity equation), and Mj = (of), 
where (•) = tr (•)p 00 . We have checked that the fixed point problem for Liouvillian 
flow (IJ5) does not admit a closed form solution of a matrix product operator form of 



small finite rank as in the weakly-symmetric case [31] , thus we had to resort to numerical 
simulations. Of course, we optimized our calculations by restricting the Liouvillean to 
diagonal blocks £|{ + i j+ i} (for p^) and £|{_i_i} (for p 2 ^). For n = 4, 6,8 and 10 sites 

I The dimension of the convex set of NESS is by one smaller than the number of distinct symmetry 
labelled fixed points 2(N — 1) due to the trace constraint tip = 1. 
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Figure 1: We plot negative current —J against spin chain length n for A = 1/2 (a), 
A = 1 (b), and A = 2 (c), and for NESSs from both symmetry subspaces, {+1,+1} 
(red/squares) and { — 1, —1} (blue/circles). Data for n < 10 are numerically exact, while 
for n > 12 the statistical/Trotter error is always smaller than the size of the symbols. 
In panel (d) we plot the case A = 2 in the log-log scale, suggesting the fast decay of 
J(n). The two eye-guiding dashed lines have slopes —2.68 (red) and —1.44 (blue). 



we used exact numerical diagonalization, while for n = 12, 14 and 16 we used a wave- 
function Monte-Carlo approach called the method of quantum trajectories, as outlined 
in Appendix A of Ref. [32] (see also Ref. [H]). An efficient Trotter expansion of the 



propagator e ' with complex coefficients [32] has been employed. The stochastic time- 
dependent Shrodinger equation has been simulated until the current was equal on all 
bonds, and the statistical error was small, both within the accuracy better than 1%. 

We consider three characteristic values of the anisotropy parameter: A = 1/2, 
A = 1, and A = 2, where the linear-response-transport in the grand-canonical ensemble 
(or in an open chain without the microcanonical constraint) has been found to be, 
respectively, ballistic [31], anomalous (super-diffusive/sub-ballistic) [35], and diffusive 
[THK22]- On the other hand, linear-response transport in the micro-canonical ensemble 
in the latter case (A = 2) has been suggested to be sub-diffusive (insulating) [37]. We 
chose the bath-coupling parameters as T = 1, /i = 0.2, which put our setup in a near- 
linear-response regime, however, we also tried other values of bath parameters and the 
results did not change qualitatively. By considering a large driving parameter close to 
the maximum /i« 1, we have found negative differential conductance for both symmetry 
subspaces (for n = 8), consistent with the behavior found in the unconstrained model 
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Figure 2: The magnetization profiles for various values of A = 1/2, 1, 2 in the {+1, +1} 
subspace, (a,b,c), respectively, and comparing different chain sizes (see the legends inside 
plots), and similar for the {—1, —1} subspace in the lower row of plots (g,h,i). In the 
midle row of plots (d,e,f) we compare, again for A = 1/2, 1, 2, and for the largest system 
size attainable, the magnetization profiles from the two symmetry sectors (red/blue). 



The spin chain with four sites is particularly interesting as the subspace { — 1,-1} 
is one-dimensional, i.e, the NESS in this subspace has to be pure state, i.e. a dark state. 
Being a pure state it cannot support current (J = 0) and has a vanishing magnetization 
on all 4 sites. Explicitly, the dark state then is of the form 

KU = 4| (loiio) -liooi}) (io) 

and does not depend on any system's parameters, A,r,/z. We did not find any dark 
states of our flow (IJ5) for n > 6. 

In Fig. [T] we display the behavior of the steady state current J versus the chain 
length n for the two distinct steady states poo and p^ in the spaces {+1, +1} and 
{—1,-1}. We find a general trend of convergence of the two curves which suggest 
that the current might not depend on the symmetry sector in the thermodynamic limit 
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Figure 3: The scaling of the spectal gap R with the system size n in the log-log scale, for 
the three values of anisotropy A = 1/2 (a), A = 1 (b), and A = 2 (c), and for the two 
symmetry subspaces, {+1,+1} (red/squares) and { — 1,-1} (blue/circles). The dashed 
lines indicate the slope —3. 



and would therefore be unique. The same applies to magnetization profiles displayed 
in Fig. |2j This implies a kind of ergodicity and suggests physical irrelevance of the 
symmetry S in the thermodynamic limit. Furthermore, in the regime A = 1/2 the 
current J(n) tends to converge to a constant suggesting a ballistic behavior, for A = 1 
J(n) seems to slowly (perhaps slightly super-diffusively) decrease with n, while for 
A = 2, J{n) decreases fast. As far as our numerics can suggest it seems that in the 
regime A = 2, the decrease of J(n) is definitely much faster than 1/n which is compatible 
with the suggested insulating behavior [37] . In Fig. [2] we make a detailed analysis of 
the magnetization profiles in both symmetry sectors, and for different system sizes. We 
find, consistently, that the profiles in the regimes A = 1/2, A = 1, and A = 2, display 
ballistic, (slightly super-)diffusive, and sub-diffusive (insulating) behavior, respectively. 
In is perhaps interesting to remark that in the isotropic case A = 1, for the largest 
system size that we could simulate, the profile seems to be almost perfectly linear which 
might also be compatible with thermodynamically normal diffusive behavior in this 
regime (perhaps related to [38J). 

In Fig. [3] we focus on the dependence of the relaxation rates of p a (t) to the 
respective steady state p^ on the chain length n. The relaxation rate is determined 
as the Liouvillean spectral gap R = minA J ^o( — Re Xj) where Xj denote the eigenvalues 
of the Liouvillian superoperator t in the respective symmetry subspace ({+1,-1-1} or 
{—1,-1}). We find consistently that in the {+1,+1} sector the spectral gap always 
quickly decays with n, perhaps faster than n~ 3 , which is the gap-scaling derived 
analytically or numerically suggested for some integrable spin chains with boundary 
driving [13, 35J. On the other hand, we are unable to conclude anything meaningful on 
the scaling of the gap in the { — 1, —1} symmetry sector, where the behavior of R(n) 
may not even be monotonic. To complement the spectral information we plot in Fig. [4] 
the set of all Liouvillian eigenvalues Xj which lie sufficiently close to the imaginary line 
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Figure 4: The full Liouvillian spectra {Xj} for the case of A = 2 and for three different 
chain sizes n — 6 (a), n — 8 (b), and n — 10 (c), and for the two symmetry subspaces, 
{+1,+1} (red/squares) and { — 1,-1} (blue/circles). We plot just a small part of 
the complex plane close to the imaginary line (i.e. plotting just the smallest rates 
corresponding to longest-lived modes). 
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Figure 5: The cumulative distribution W(r) of the the case A = 2, n = 10, i.e. the 
probability that the randomly chosen decay rate —Re Xj is larger than — r, for the two 
symmetry subspaces, {+1, +1} (red) and { — 1,-1} (blue). Note the occurrence of two 
pseudo-gaps (strips in complex plane with very few - rarely distributed decay modes), 
particularly visible in the data for { — 1,-1}. 



(i.e. with sufficiently small damping rates — ReAj). Further, we considered the density 
distribution of Liouvillian spectrum projected onto the real line, in fact a cumulative 
distribution W{r) = J2\ $( r — R e -\?)/^A 1 giving the probability that a randomly 
picked damping rate is larger than — r. Here 8(r) designates a Heaviside step function. 
The plot of W(r) for the two symmetry sectors shown in Fig. [5] in the easy-axis regime 
A = 2 suggests existence of two pseudo-gaps, particularly clearly for the { — 1,-1} 
symmetry sector. 
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3. Discussion and conclusion 

We have presented a small but potentially very useful observation on the symmetry 
reduction of the Linbdblad master equations describing Markovian open quantum 
systems. We classified quantum Liouvillian symmetries as weak or strong, depending, 
respectively, on whether only the generator of master equation as a whole commutes 
with the symmetry operation, or whether each of the Lindblad (jump) operators, and 
the Hamiltonian, commute with the symmetry individually. We have shown that 
only existence of a strong symmetry implies nontrivial symmetry reductions inside 
the space of steady states (fixed points of the Markovian semi-group). We have 
also provided a non-trivial example of our construction in the context of quantum 
transport problem for the strongly interacting XXZ spin 1/2 chain. Namely defining an 
open XXZ chain with a micro-canonical constraint, we have shown that it exhibits a 
strong parity-type symmetry and outlined numerical computation of distinct steady 
states and the corresponding physical observables. Our numerical results seem to 
suggest that microcanonically constrained open XXZ chain exhibits a sub-diffusive 
(thermodynamically insulating) behavior in the gapped Ising-like (easy axis) phase. 

We note that both, weak and strong symmetries can provide practical advantage 
in description of quantum Liouvillian flows, as they both allow to block-diagonalize 
and hence reduce the dimension of the Liouvillian. One can also have a combination 
of both, namely independent generators of weak and strong symmetries, or non-abelian 
symmetries of either type. In particular, the study of non-abelian Liouvillian symmetries 
should be an important line of future research. 
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4. Appendix: Symmetry reductions of quantum Liouvillian dynamics 

Let us consider two cases of symmetric Lindblad master equations with two different 
kinds of dissipators of the flow (jl]) : 

(i) For the first kind, which we will call a weak symmetry, we assume that there exist 
a unitary operator S over TL, such that 

£{SxS j ) = S(Cx)S j (11) 

for any x G £>("H). 
(ii) For the second kind, which we will call a strong symmetry, we assume that a unitary 
operator S over % commutes with each member of the set {H, L m , m = 1,2,...} 
individually 

[S,H} = 0, [S,L m } = 0,m = l,2,... (12) 
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Clearly (ii) implies (i), but not vice versa. 

We continue by fixing some notation. Let s a = e l9a , a — 1, 2 . . .715 denote distinct 
eigenvalues of S, and let 7i a be the corresponding mutually orthogonal eigenspaces, so 
that we have a complete symmetry decomposition of the Hilbert space 

I'S 

H = 0W„. (13) 

a=l 

The adjoint representation of S on Hilbert space B(H) shall be denoted as S, 

S{x) = SxSl (14) 

The spectrum of a super-operator S is formed of all possible products s a S/3 = e l( - ea ~ e0 \ 
since S (\if;) ((f)\) = s a sp\il>)((f)\ for any \ip) <G 7i a , \(p) e %p. Let s' u , v — 1, 2, . . .n s denote 
the distinct eigenvalues of S, denoted such that s[ = 1, and s' v ^ 1 for all v ^ 1. Note 
that the following bounds always hold 

n s < n§ < n s (n s - 1). (15) 

The lower bound ns = ns is reached when {s a } are roots of unity, i.e. in the case of Z ns 
symmetry obeying S ns = 1. The upper bound n s = ns(ns — 1) is reached when there 
is no arithmetic structure in the spectrum {s a } which may happen when S generates a 
continuous group of symmetry transformations. Let 

n s 

B(U) = @B U (16) 

be the corresponding symmetry decomposition of the operator space B(H), where the 
5-eigenspaces B v are mutually orthogonal in the Hilbert-Schmidt sense. 



Clearly, having a weak symmetry (11) is equivalent to having commuting 
superoperators 

is = SC (17) 

which means that the Liouvillian is block diagonalized with respect to the decomposition 



(16), namely 

£B V C B v , (18) 

i.e. all spaces B v are invariant w.r.t. the flow (JT1). Note that only B\ contains operators 
with non-vanishing trace, thus it should contain the fixed point p^ e B\ due to the trace 
preservation. We remark that under the conditions of the Evans theorem [22] NESS can 
be unique despite the existence of a weak symmetry, as shall be the case in the example 
outlined in the next section. 

Let us now assume the existence of a strong symmetry. In this case we can prove 
the following useful result: 



Theorem 4.1 Existence of a unitary operartor S with the property (12) implies: 
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(i) Liouvillian C can be block-decomposed into rig invariant subspaces 

where B afi = {\ip) {(f>\; \ij)) G U a ,\4>) e %} f or a > P = 1,2, ...715. 
(ii) We have at least ns distinct fixed points (NESSs), which we can label by s a , namely 
each diagonal operator space contains at least one fixed point, 

Poo G l5a,a 

for a = 1,2, ...ns. 

Proof Let us write the super-operators of left and right multiplications with the 
symmetry operator S^(x) = Sx, Sr(x) = xS\ for any x G B(7i). Since, clearly, 
[S'l, Sr] = 0, the spaces B a ^ can be considered as the joint eigenspaces of both Sl and 
S R . Namely S L (|V>)<0|) = L\iJ>)(<f>\, 3 K (\if>) (<f>\) = sp^ttl for an Y 1^) e K a , \<j>) G Up. 



Eq. (12) implies commutativity [.£, S'l] = 0, [C, Sr] = 0, implying that any B a ^ is also 
invariant subspace of C proving point (i) of the theorem. 

To prove (ii) we note again that operators with non-vanishing trace can only be 
contained in diagonal spaces B a , a , due to mutual orthogonality of T-L a . Hence starting 
from some p a (0) G £> Q)Q with trp a (0) = 1 (say p(0) = P a / tr P a where P a is orthogonal 
projector to the subspace Ha), the flow (fTl) yields a fixed point p^ in B a>a for any 
a = 1, . . . , ns- Again according to Evans theorem [22] this fixed point can be unique, or 
further degenerate. In any case, we have at least n s distinct steady states, which may 
be labelled by the symmetry eigenvalue s a . QED 

Clearly, the direct sum of all invariant subspaces give the entire operator space 

n s n s 

w = 00*w> ( 19 ) 

a=l /3=1 



and the quotient invariant spaces ( 16 ) are given by partial direct sums 

Sa sp=s l v 

B v = B a3 . (20) 

We remark that (non-positive) initial conditions p(0) supported in the non-diagonal 
invariant spaces B a< p, a ^ (3, which have vanishing trace, could in principle provide extra 
degeneracy of the null space of C, i.e. an extra dimension of the convex set of NESSs, or 
simply complete the description of the relaxation process for a general initial condition 
p(0). In other words, we do not see any argument which would forbid the spectra of 
block-operators C\s a to contain even for a ^ (3, although this should be - if possible 
at all - an exceptional (non-generic) situation. 
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